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The  first  Bom  Approximation  is  used  to  j: alculate  the  differential 
cross  section  for  electron  detachment  from  H~  ions  colliding  with 
hydrogen  atoms.  The  differential  and  total  cryss  sections  are  evalu¬ 
ated  for  the  resultant  hydrogen  atom  in  the  Is  or  2s  state  and  the 
target  hydrogen  atom  in  any  ecited  state  up  to  n  =  8.  The  scattering 
matrix  includes  the  exchange  of  particles  interaction,  which  has  not 
been  included  in  earlier  work.  The  results  are  presented  and  compared 
with  earlier  work  where  possible. 
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DIFFERENTIAL  CROSS  SECTION  FOR  ELECTRON 
DETACHMENT  FROM  NEGATIVE  HYDROGEN  IONS 
COLLIDING  FROM  NEUTRAL  HYDROGEN  ATOMS 

I.  Introduction 

Background  and  Previous  Work 

The  recent  interest  in  neutral  particle  beams  has  made  it  necessary 
for  investigators  to  better  understand  the  theoretical  and  experimental 
processes  in  the  production  of  such  beams.  There  are  two  basic  pro¬ 
cesses  that  can  be  used  to  create  a  beam  of  neutral  particles.  The 
first  method  is  to  pass  a  beam  of  positively  charged  ions  through  a 
gas  pick-up  cell,  where  the  ions  pick-up  an  electron.  The  resulting 
beam  consists  of  neutral  particles.  The  second  process  is  to  pass  a 
beam  of  negatively  charged  ions  through  a  gas  stripping  cell,  stripping 
the  extra  electron  from  the  particles  in  the  beam.  The  stripping  or 
pick-up  reaction  has  deflected  the  outgoing  beam  of  particles  because 
of  a  change  of  energy  and  hence  there  is  a  momentum  transfer.  The 
stripping  reaction  would  deflect  the  beam  less  because  of  the  smaller 
change  in  energy  of  the  atomic  systems  as  compared  to  the  pick-up  reac¬ 
tion.  Using  the  stripping  reaction  a  beam  of  neutral  hydrogen  could  be 
produced  by  stripping  an  electron  from  a  negatively  charged  hydrogen 
ion  (H_)  in  a  collision  with  hydrogen  or  a  heavier  atomic  system. 

There  are  many  theoretical  questions  concerning  the  process  of 
electron  detachment.  What  is  the  differential  cross  section  for  differ¬ 
ent  energies?  What  effect  do  better  wavefunctions  have  on  the  theoreti¬ 
cal  differential  cross  section?  How  important  is  exchange?  Some  of 
these  questions  have  been  answered  in  the  literature.  The  basic 


1 


theoretical  concepts  for  electron  detachment  were  presented  by  Bates  and 
Griffing  (Ref  2  and  3),  however  their  work  dealt  with  finding  the  total 
cross  section.  Sida  (Ref  17)  calculated  the  total  cross  section  for 
electron  detachment  of  H“  ions  with  helium  atoms.  The  process  of  H~ 
ions  colliding  with  H  atoms  was  done  by  McDowell  and  Peach  (Ref  15). 

In  both  cases  simple  two  parameter  wave  functions  were  used  and  no  ex¬ 
change  of  particles  was  discussed. 

Recently  there  has  been  a  greater  interest  in  obtaining  the  dif¬ 
ferential  cross  section  because  it  shows  explicitly  the  deflection  of 
the  beam  from  its  original  direction.  Lee  and  Chen  (Ref  13)  have  cal¬ 
culated  the  differential  cross  section  for  electron  detachment  from  H” 
ions  in  collisions  with  H  atoms  using  a  sum  rule  technique  and  summing 
over  all  the  excited  states  of  the  resultant  H  atom.  Genoni  and  Wright 
(Ref  8)  have  also  calculated  the  differential  cross  section  for  this 
process.  Their  calculations  were  for  the  target  H  atom  in  final  state 
N  =  1  to  8  or  in  a  continuous  state  and  the  resultant  H  atom  in  the  Is 
or  2s  state.  Exchange  of  particles  was  neglected  in  both  cases. 

Purpose  and  Scope 

The  purpose  of  this  thesis  is  to  calculate  the  differential  cross 
section  for  electron  detachment  from  H”  ions  colliding  with  H  atoms. 

The  differential  cross  section  is  evaluated  using  the  first  term  of  the 
Bom  Series.  The  H-  ion  is  initially  in  its  ground  state  configuration 
and  has  a  kinetic  energy  of  1  to  500  MeV  in  the  laboratory  frame,  while 
the  target  hydrogen  atom  is  in  its  ground  state  configuration  and  sta¬ 
tionary  in  the  laboratory  frame.  After  the  collision  the  resultant  H 
atom  is  in  the  Is  or  2s  configuration  while  the  H  target  atom  is  in  the 
N  =  1  to  8  state.  The  differential  cross  section  and  total  cross  section 
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is  evaluated  explicitly  for  energies  of  2  MeV  and  200  MeV  for  the  H”  ion. 
Units 

All  quantities  in  this  study  are  presented  in  atomic  units,  where 

■h  =  m  =  e  =  1 

and  are  dimensionless.  The  speed  of  light  is  the  inverse  of  the  fine 

structure  constant,  a.  Length  and  momenta  are  also  dimensionless  while 

energy  is  in  units  of  hartrees,  where  1  hartree  =  27.21  eV.  All  cross 

2 

sections  are  presented  in  units  of  7raQ  where  aQ  is  a  Bohr  radius. 

General  approach 

The  general  approach  in  the  evaluation  of  the  differential  and 
total  cross  section  follows  the  pattern  set  in  solving  any  quantum 
mechanical  problem.  The  general  kinematics  are  discussed  with  a  spe¬ 
cial  emphasis  on  the  relativistic  treatment  of  the  scattering  problem. 

The  total  Hamiltonian  is  formulated  for  the  system  of  particles  and 
from  the  total  Hamiltonian  the  scattering  potential,  initial  and  final 
wavefunction  for  the  particles  are  found.  The  scattering  matrix  can 
be  evaluated  using  the  scattering  potential  and  wavefunctions .  Using 
this  scattering  matrix,  the  differential  and  total  cross  section  are 
obtained  using  numerical  integration  techniques. 
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II.  Formulation  of  Scattering  Problem 


In  this  chapter  the  basic  scattering  problem  is  discussed  and  ex¬ 
plained.  An  overview  of  the  scattering  is  presented  with  special  detail 
given  to  the  energy  and  momentum  considerations.  The  differential 
cross  section  for  the  final  three  body  problem  is  derived  using  Fermis' 
Golden  Rule.  Finally  a  discussion  of  the  Born  Series  is  given  and  it  is 
shown  that  only  the  first  term  of  the  Born  Series  needs  to  be  retained 
in  the  calculation  of  the  scattering  matrix. 


Overview  of  the  Scattering  Problem 

The  scattering  problem  in  general  is  that  of  an  incoming  H~  ion, 
with  energy  T,  traveling  toward  a  stationary  H  target  atom.  After  the 
collision  there  are  three  particles  in  the  final  state,  a  free  electron 
and  two  hydrogen  atoms.  The  resultant  H  atom  (H-  stripped  of  one 
electron)  has  been  deflected  through  some  scattering  angle  0 .  The 
differential  cross  section  associated  with  the  scattering  angle  0  is  to 
be  calculated. 

The  differential  cross  section  can  be  evaluated  in  the  laboratory 
frame,  but  it  is  easier  to  perform  the  calculations  in  the  center  of 
mass  frame  and  then  transform  back  to  the  laboratory  frame.  Therefore 


in  the  center  of  mass  frame  the  H”  ion  has  momentum 

k  =  0 _ k 

-  M_  +  M  i 


(1) 


where  MQ  is  the  mass  of  the  neutral  H  atom,  M_  is  the  mass  of  the  H" 
ion  and  is  the  initial  momentum  of  the  H~  ion  in  the  laboratory 
frame.  When  relativistic  effects  are  included  k^  is  given  by 


T2  +  2TM_c2  v± 


(2) 


Mm* - maim* 
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where  T  is  the  incoming  energy  and  is  the  unit  vector  of  the  initial 
velocity.  This  equation  is  valid  for  all  incoming  energies  that  need  to 
be  considered.  In  the  center  of  mass  frame  the  H  target  atom  has 


momentum 


*o  =  -K- 


The  total  energy  of  the  initial  system,  in  the  center  of  mass  frame,  is 

k2  k2 

Ea  =  2FT  +  2M^  +  ( 1s  ]  +  eh  ( 1s)  ^ 

where  E^-  (Is  )  is  the  energy  bound  in  the  H  ion  system  and  EH  (Is) 
is  the  energy  bound  in  the  H  target  atom  system.  E^  (Is)  is  given  to 


Ey  (Is)  =  -  0.5  a.u.  (5) 

and  from  Caveline,  et  al.  (Ref  6:625), 

E^  (Is2)  =  -  0.527751  a.u.  (6) 

Using  Equation  3,  Equation  4  can  be  rewritten  as 

k  2  » 

E0  =  —  +  E„  (Is)  +  E„-  (is*)  (7) 

a  2y  "  ^ 

-o 

where  u  is  the  reduced  mass  of  the  H~  ion  and  H  target  atom. 

-o 

The  final  configuration  consists  of  three  particles,  a  free  elec¬ 
tron  with  momentum  k  ,  a  recoiling  target  atom  with  momentum  1c  and  a  re- 

C 

suiting  H  atari  with  momentum  k^.  From  the  conservation  of  momentum 
argument,  the  momentum  of  the  H  target  atom  is 

k  =  -  (kjj  +  kg)  (8) 

The  total  energy  of  the  final  system  is  then 
1<2  ^  k  2 

%  *  a-  -  3T +  sr  *  %  (nl)  *  %  ,N)  *  hr  (9) 

o  o  e 

where  E^  (nl)  is  the  energy  bound  in  the  resultant  H  atom  when  nl  =  Is 


or  2s,  (N)  is  the  energy  bound  in  the  target  atom  when  its  final 

state  is  N  =  1  to  8  and  1^-  is  the  ionization  energy  of  the  H“  ion, 


from  Ref  6, 


IH-  =  0.76  eV  =  2.793  x  10-£  a.u. 


Using  Equation  8,  Equation  9  can  be  rewritten  as 

\r  2  2  2k 

Eb  =  ar +  TT +  -ar2  +  %  (N)  +  ^  {nl)  +  V  l11) 

o  Hoe  o 

where  y  is  the  reduced  mass  of  a  hydrogen  atom  and  an  electron, 
oe 

The  conservation  of  energy  requires  that 


+  Ejj-  (is*)  +  EH  (is)  = 


ke2 


+  Ejj  (N)  +  EH  (m)  +  In- 


Denoting  _ 

k  * 

Ef  =  +  Ejj"  (IS2)  +  Ejj  (IS)  -  Ejj  (til)  -  Ejj  (  N  )  -  Ijj-  (14 

k  2 

-  ll 

=  - - 1.0  a.u.  +  .5  —o  +  ~o  a.u.  (15 

2^  n2  N2 

allows  Equation  13  to  be  written  as 

Kh2  ... 

f  =  Mo  +  2yeo  + 

The  term  on  the  right  represents  all  the  known  quantities  while  the 
terms  on  the  left  represents  the  unknowns  kjj  and  kg. 

Formulation  of  Differential  Cross  Section  With  Three  Final  Particles 
The  scattering  of  one  particle  with  another  can  be  treated  as  a 
time-dependent  perturbation  where  the  potential  is  "turned  on"  when 
the  incident  particle  is  in  the  vicinity  of  the  target  particle,  and 
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then  the  potential  is  "turned  off"  as  the  particles  leave  the  vicinity 
of  each  other.  The  scattering  process  can  be  considered  as  a  time- 
dependent  transition  from  an  initial  wave  state,  | i>  to  a  final  wave 
state,  |f>.  For  the  problem  under  consideration,  the  initial  wave 
state  consists  of  the  wave  functions  for  the  H  target  atom  and  the  H” 
ion.  The  final  wave  function  consists  of  the  wave  functions  for  the  H 
target  atom,  the  resultant  H  atom  and  the  free  electron.  The  transition 
rate  per  incident  particle  is  given  by  Fermis '  Golden  Rule 

R  =  |<f|v|i>|2  p(Ef)  (17) 

where  V  is  the  scattering  potential  and  p(E^)  is  the  density  of  final 
states.  The  integral  < f J V | i>  are  the  scattering  matrix  elements  and  will 
be  denoted  as  T^.  The  density  of  final  states  can  be  determined  for 
three  particles  using  a  phase  space  argument. 

If  there  are  three  final  particles  of  mass  n^,  momentum  k^,  and 
obey  the  constraint 

Tc.j  +  Tc2  +  lc3  -  0  ( 18) 

then  the  system  has  6  degrees  of  freedom.  Assuming  Ic^  and  are  the 

independent  momenta,  the  number  of  states  in  an  interval  of  momenta 

die.  die-  is  , 

lP  dk.  dk5 

dn  =  - '-z— -  (19) 

(2tt)° 

3  3 

where  L  is  a  volume  normalization  constant,  usually  equal  to  (2tt)  . 

Planck's  constant,  h  does  not  appear  in  Equation  19  because  "h  =  1  in 

atomic  units.  The  density  of  states  is  defined  as  the  number  of  states 

per  energy,  or 

»<¥  =  §rt  t20) 

where  E^.  is  the  total  energy  of  the  final  state  configuration,  and  is 
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given  by 


k  2  k  2  k  2 

-  _  1  ^  K3 

f  ~  2m1  +  2m2  +  2ni2 

Conservation  of  energy  requires  5 

t  3  ki 

7s  <Er  -  J,  ar)<JEf  *  1 

Multiplying  Equation  22  by  dn  and  integrating  gives 

*  =  £ dn -4^  l  di?26(Ef -Is: 


(2tt)  AE 


1231 


This  integral  represents  the  total  number  of  states  in  the  energy  range 
AEf  when  all  possible  combinations  of  dl<2  and  dl<2  are  taken.  In  the 


limit  as 


limit  An  _  dn 
AE^->o  AE^  “  dEj. 


implies  that 


<=<¥  =  *2  SIEf  '  j,  ST’ 

12,)6  ^  1=1  1 


which  can  be  written  as 

r  6  r 


3  k. 


p(Ef)  =  (^r)  /k,2  dk1  dn,  k22  dk2  dn2  5(Ef  -  J  ^-)  (26) 

It  is  now  possible  to  integrate  over  dk2  space  using  the  properties  of 
the  delta  function.  The  argument  of  the  delta  function  can  be  re¬ 
written  using  Equation  18,  where  ?  ? 

k,  k2  2lc,  ’1<2 

f(k1fk5)  =  Er  -  -5-! - *= - 5^; — —  (27) 

12  f  2u13  ^23  ^3 

where  ^  is  the  reduced  mass  of  particles  i  and  j.  Using  the  identify 

S(f(k1-k2»  =  J,  |i"(kX>|2  ' -  (28> 


k2  =  k2 


where  there  are  n  roots  of  f ( k, ,k2)  and  k2  is  a  root,  so  that 

f(k.,,k5a)  =  0 


8 


The  roots  of  f(k.j,k2a)  are  found  to  be 
k^  cos  \J) 


"5  ' 


“3 


Ik2  cos24>  2eI k  2 

+  / - 2 - + - 

1  m3  y23  y23y13 


23 


(30) 


where  \jj  is  the  angle  between  ^  and  k^.  The  negative  root  can  be  ig¬ 
nored  because  it  is  impossible  for  the  magnitude  of  a  real  vector 
quantity  to  be  negative.  Therefore  there  is  only  one  root  for  f(k1,k2), 


m- 


v 


2  2  y 

k^  cos  lk^  cos  2Ef  k.  ”* 

+  ./ - - - + - 

y23  y23y13 


“3  f  m3 

The  derivative  of  f(k^,k2)  is 


23 


(3i : 


f’(krk2)  = 


df(krk2) 


dk2 

k2  k^  cos 
y23  m3 


therefore 


f'(krk2) 


k2  =  k2+  ^ 


k^2  cos2  2Ef  k2 
2 - + - 


(32) 

(33) 

(34) 


"3 


y23  y23y13 


Performing  the  integration  of  k2,  using  Equation  28  gives  the  density 


of  states  to  be 

k2  dk1  dS21  dfi0  yiL 
p(Ef)  =  — - - —  1  2  2J 


k^  cos  \(>  ^^1^7 


m3 


m-j 


2  cos2^  2Ef  k2 


±23  y22y.1.3. 


*1  ' 

i 


(35) 


*3  y23  y23y13 

The  transition  rate  per  unit  solid  angle  due  to  a  particle  flux  of 
intensity 

V, 

(36) 


-3 


is 


R  T  da 
dfl  ‘  x  dft 


(37) 


The  differential  cross  section  is  then  given  by 
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where  p(Ef)  is  given  in  Equation  35. 

For  the  specific  problem  to  be  evaluated  the  following  substitu¬ 
tions  are  made 

*1  =  ke 

k2  =  1<H 
“2  =  m3  =  m0 
m1  =  me  =  1 

and  Ef  is  given  by  Equation  15. 


The  Born  Approximation 

The  scattering  matrix  elements,  T_^,  associated  with  the  differ¬ 
ential  cross  section  depends  upon  the  initial  and  final  wave  functions 
and  also  the  scattering  potential.  If  the  scattering  potential  is  dif¬ 
ferent  before  and  after  the  scattering,  the  potential  can  be  written 
as  a  series,  where 


V  =  V  +  V  G  +  +  V  G  *  VK  G+  Va  +  - 

a  a  o  b  aoboa 


(39) 


where  Va  is  the  scattering  potential  before  the  collision,  Vb  is  the 
scattering  potential  after  the  collision  and 


p  +  1 

Go  '  E_  +  in  -  H_ 
a  o 

where  GQ+  is  a  form  of  a  Green’s  function.  This  expansion  of  the 
potential  gives  rise  to  the  Born  series  which  is  given  by 
<f|V|i>  =  <f|Va|i>  +  <f|Va  G0+  Vb|i>  +  - 

where  if  all  but  the  first  term  of  the  expansion  are  ignored,  or 

<f|V(i>  *  <f|Vali> 


(40) 


(41) 


(42) 


is  called  the  first  Born  Approximation.  If  the  second  term  is  included 
in  the  expansion  it  is  then  called  the  second  Born  Approximation,  etc. 

The  first  Born  Approximation  is  valid  as  long  as  the  incoming 
energy  of  the  particle  is  much  greater  than  the  binding  energy.  In 

L  g 

this  study  WQ  =  0.5  har trees,  while  «  10  -  10  har trees,  therefore 
the  first  Born  Approximation  can  be  used. 

At  high  energies  the  second  term  of  the  can  be  approximated  as 


T'2’  *  r-  <flva2l1> 


(43) 


An  estimate  of  the  correction  supplied  by  the  second  term  is  found  by 
looking  at  the  ratio  of  to  ,  or 


and  in  all  probability 


Then 


t(2) 

1 

<f|V|i> 

7HT 

=  Ei 

<f|Va|i> 

io6, 

t(2) 

1 

<f|Va2|i> 

^TF 

<f|Va|i> 

<f|V|i>| 

<f|Va|i>l 


~  10 


|T{2)|  -  10"V(1)| 


(44) 


(45) 


(46) 


(47) 


and  the  first  term  of  the  Bom  Series  can  be  used  to  calculate  the 
differential  cross  section,  accurately. 
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III.  The  Hamiltonian,  Scattering  Potential, 
and  Particle  Wave  Functions 


The  motion  and  subsequent  wave  functions  for  a  system  of  particles 
can  be  found  by  solving  the  quantum  mechanical  Hamiltonian  of  the 
system.  This  chapter  develops  the  Hamiltonian  of  the  total  system 
and  from  it  the  scattering  potential,  and  initial  and  final  wave  func¬ 
tions  are  obtained.  This  chapter  concludes  with  a  discussion  of  the 
antisymmeterization  of  the  final  wave  function. 


The  Total  Hamiltonian 

The  motion  of  the  particles  in  the  center  of  mass  frame  is  des¬ 
cribed  completely  by  the  total  Hamiltonian  of  the  system.  In  the  case 
of  a  H"  ion  colliding  with  a  H  atom,  the  Hamiltonian  is  given  by 


1 


'1 


P2 


H  =  -3r  +  sr  +  sr  +  3r  +  sr" 

p  p 


_ 1 _ 1 - 2_  +  _L_  + 

(R.J  -  X,|  |^2  “  ^2 1  1^2  ”  ^3 1  1^2  ”  ^3 1 

I*!  -s2|  pr,  -  x2i +  -x3|  |^-X2| 

- L_ - 3 -  (48) 

l*,  -X3I  1^2  “  *1 ' 

where  the  first  two  terms  represent  the  kinetic  energy  of  the  two  pro¬ 
tons,  the  next  three  terms  represent  the  kinetic  energy  of  the  electrons 
and  the  last  ten  terms  represents  the  electrostatic  potential  between 
the  particles.  The  position  vectors  in  the  electrostatic  potential 
terms  are  depicted  in  Figure  1 ,  where  is  the  position  vector  of  the 
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proton  and  is  the  position  vector  of  the  electron.  The  total 
Hamiltonian  can  now  be  used  to  find  the  particle  Hamiltonians  and  the 
scattering  potential  before  and  after  the  collision. 

When  the  particles  are  separated  and  do  not  interact,  the  initial 
configuration  consists  of  a  H~  ion  and  a  H  atom.  The  Hamiltonian  for 
the  H  atom  is  given  simply  by 


1 

I*,  "  *-,1 


(49) 


assuming  electron  1  and  proton  1  are  in  the  H  target  atom  initially. 
The  H”  ion's  Hamiltonian  is  the  three-body  Hamiltonian 


2me  2me  IX3-X3I 


1 _ 1 

|I2  -  x2j  |R2  -  X3| 


(50) 


assuming  proton  2  and  electron  2  and  3  are  in  the  H“  ion.  In  reality 
the  designation  of  each  particle  is  for  accounting  purposes  only,  and 
comes  into  play  when  there  is  an  exchange  of  particles. 

The  scattering  potential,  V  ,  can  be  found  by  subtracting  the 

Cl 

particle  Hamiltonians  from  the  total  Hamiltonian  or 


Va  =  H  -  %  -  V 

=  — L - +  — — 1— —  +  — — 1— — 

"  ^2 1  lxi-x2l  lxi  -  X3I 

_ 1 _ 3 _ 1 

|In-X2S  |R,  -  X3| 


(51) 


(52) 


To  simplify  notation,  the  following  convention  will  be  used 


(53) 
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(54) 


V  i  =  - - - 

J  \\  -  y 

v  - - 3 -  (55) 

P  IB,  -F2I 

The  scattering  potential  is  now  written  as 

va  =  Vp  *  V,2  *  V13  *  V  -  V  -  V12  1561 

After  the  collision,  when  the  particles  are  separated  and  do  not 
interact,  there  are  two  H  atoms  and  one  free  electron.  The  Hamiltonian 
for  the  H  target  atom  is  o  o 

•  V  pi  a 

*h  =  sr +  sr "  vi  (57) 

P  e 

tu  Ll. 

where  a  represents  the  a  proton  and  i  represents  the  i  n  electron. 

The  Hamiltonian  for  the  resultant  H  atom  is 

where  proton  B  and  electron  j  are  in  the  resultant  H  atom.  The 

Hamiltonian  for  the  free  electron  k  is 

i  pif 

Kq  =  g-  (59) 

e 

The  above  notation  was  used  to  account  for  exchange  of  the  particles 
during  the  collision.  In  the  above  cases  if  a  =  1,  B  =  2  or  if  B  =  1, 
a  =  2,  and  i,  j,  and  k  are  the  permutations  of  1 ,  2  and  3.  All 
permutations  of  the  indices  are  presented  in  Table  1 .  The  permutation 
of  the  particles  is  discussed  in  further  detail  when  the  final  wave 
function  is  formulated. 

Once  again  the  scattering  potential  is  found  using 

vb  *  H  -  *H  -  S  *  K  '  l60) 

=  \  ♦  *  vik  *  vjk  - -  V-  vi6-  Vk6  1611 
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Table  1 

Permutations  of  Particles 

in  Final  State 

H_ 

(Target) 

H  (Resultant) 

e"  (Free) 

e~ 

P 

e~ 

P 

e~ 

1 

1 

2 

2 

3 

1 

1 

3 

2 

2 

2 

1 

1 

2 

3 

2 

1 

3 

2 

1 

3 

1 

1 

2 

2 

3 

1 

2 

2 

1 

1 

2 

2 

1 

3 

1 

2 

3 

1 

2 

2 

2 

1 

1 

3 

2 

2 

3 

1 

1 

3 

2 

1 

1 

2 

3 

2 

2 

1 

1 

Therefore  Va  i  Vb,  that  is  the  scattering  potentials  are  different 
before  and  after  the  collision.  It  can  be  shown  that  the  scattering 
matrix,  T^,  is  the  same  whether  Va  or  is  used.  Tfi,  is  given  by 

Tfl  »  <f|»a|i>  <« 

Using  Equation  51  this  can  be  written  as 

Tfl=  <f|(H  -  Kn  -  K^-)i>  (62 

and 

(Kn  +  Kn-)  |i>  =  Ea|i>  (64 

therefore  Equation  63  becomes 

Tfi  =  <f | (H  -  Ea)  i>  (65 
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(66) 


Since  H  and  E_  are  both  hermitian  operators  this  can  be  written  as 

Tfi  =  <(H  -  Ea)f|i> 

But  on  the  energy  shell  Ea  =  Eb  (Conservation  of  Energy) , 

Tfi  =  <(H  -  |i>  (67) 

and  using  Equation  60  this  is 

Tfi  =  <(H  -  -  Kg  -  Ka)f|i>  (68) 

or 

Tfi  =  <Vbf|i>  (69) 

Since  Vb  is  also  hermitian,  Equation  69  is 

Tn  =  <f|Vb|i>  (70) 


Therefore 

Tfi  =  <flvali>  = 

and  the  scattering  potential  will  be 


=  VK  +  V,-  +  V„  -  V 


'12 


13 


(7D 


(72) 


Initial  Wavefunction 

Hie  initial  particles  wave functions  are  obtained  by  obtaining 
solutions  to  the  particle  Hamiltonians.  The  initial  wavefunction  is 
the  product  of  the  initial  particle  wavefunctions .  The  initial  wave 
function  of  the  H  target  atom  is 


\l>H(1s,1)  =  Nh  exp(-r.J)exp(-iIc_*T?0)Uv  (s.pu^  (P^ 


(73) 


where  NH  is  a  normalization  constant  equal  to  (8n  )  the  first  term 
represents  the  atomic  wavefunction  for  an  electron  in  the  ground  state 
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of  a  H  atom,  the  second  term  represents  the  plane  wave  associated  with 
the  center  of  mass  of  the  H  atom,  the  third  and  fourth  terms  represent 
the  spin  functions  for  the  electron  and  proton,  respectively.  Equa¬ 
tion  73  written  in  parabolic  coordinates  is 


%( Is;  1 )  =  Nh  expl-Ju  +  q))exp(-ik*Ro)Uv  (s.,)Uy  (P^ 


(74) 


The  wave  function  for  the  H  ion  is  a  numerical  solution  to  the 
three  body  Hamiltonian  (Equation  50).  There  are  many  numerical  solutions 
in  the  literature.  For  this  study,  a  simple  two  parameter  wavefunction 
will  be  used  (Ref  6:625)  in  order  that  the  spatial  integration  to  be 
done  analytically.  The  H“  ion  is  given  by 


2 

4^-(1s  ;2,3)  =  Nh-  [exp(-ccr2  -  +  exp(-ar3  -  Br2f]  • 

exp ( ik  -R  )  [U,  (sJU  n  (s~)  -  Un (s„)U  n(s,Q-' 
-  -  S  3  -5  J  -z  k 


yv 


(75) 


fa3B3 


where  N„-  is  a  normalization  constant  equal  to  the  first  term 

f  32ir’J 

in  brackets  is  the  orbital  wave  function  for  the  H"  electrons  with 
a  =  1.03923  and  B  =  0.283221,  the  second  term  represents  the  plane  wave 
associated  with  the  center  of  mass  of  the  H"  ion,  the  third  term  repre¬ 
sents  the  antisymmetric  spin  functions  of  the  two  electrons  while  the 
fourth  term  represents  the  protons  spin  function.  The  total  initial 
wave  function  is  then 


|  i>  =  ipH(  Is:  1 )  \j^H(  Is  .2,3) 


(76) 


Final  Wavefunction 

After  the  collision  the  three  particles  have  wavefunctions  that 
correspond  to  the  solution  of  the  corresponding  particle  Hamiltonian. 
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A.L. 

The  H  target  atom  in  the  n  n  excited  state  has  a  wave  function,  in 
parabolic  coordinates,  given  by 

NH  1 

^(N;i,a)  =  -g  exp(-  g  (5  +  n)/N)F{-n1 , 1  ,£/N)  • 

F(-n2,1,ry/N)exp(-i(kH  +  kg)  ’^q)Uv.  (si} 


\{h> 


(77) 


.2,-1 


where  =  (8u  )  ,  F(-n^ , 1 ,C/N) ,  and  F(-n2, 1 ,n/N)  are  confluent  hyper¬ 

geometric  functions  and  N  =  n^  +  n2  +  1 .  The  resultant  H  atom  has  a 
wave  function  of  the  form 


IfcdsjJ.B)  =  Nr  exp(-rj)exp(ikH*IH)Uv,(sj)Uyf(PB)  (78) 

»  p  -1 

where  NR  =  (8u  )  2 ,  if  the  final  state  is  nl  =  Is,  or 

4»R(1s;j»B)  =  Nr(2  -  r j ) exp ( -r  j/2 ) exp ( ik^ •  RR )  • 

V’j’V1’61  (79) 

where  NR  =  (256tt2)-2  if  the  final  state  is  nl  =  2s.  The  wave  function 
for  the  free  electron  k  is  simply 


*e(k)  =  (2* )'3/2  exp(ike.re)*Ur’(sk) 


(80) 


Because  the  particles  could  have  exchanged  positions  it  is 
necessary  to  antisymmeterize  the  final  wave  f’mction.  Therefore  the 
total  final  wave  function  is 


tyf  =  Si? 


(^1 » V-|  J^i  >  Mi  ) 

(rV  v1JP2,y2  ^ 

( r2  *  v2  i  >  ) 

(r2’  v2>  P2,y2  } 

{r2,v^} 

( r3 ’ v3  ;P1,y1  5 

( r3 »  v3 » :  ^2 '  ^2  ^ 

^r3»v3> 
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/T2 


(81) 


( **i  t  V-j 

( > \>i  jP-j » y-j  ) 

( r*2 * v2 

,P2’U2  ^ 

(r2,V2  ;pvyi 

(r2>v2* 

(r3’v3 

5  ^2  *  ^2  ^ 

•1 

( r3 * V3  >Pi »y-|  ) 

{r*3 » \>2 } 

t  f 

where  (rj_>  vi5pa»^a  )  is  the  wavefunction  of  the  H  target  atom  with 

t  i  » 

electron  i  and  proton  a,  (r^,  v^;  P^,  )  is  the  wavefunction  of  the 

resultant  H  atom  with  electron  j  and  proton  3,  and  {r^, 

wavefunction  for  the  free  electron  k.  Using  the  scattering  potential, 

V,  the  initial  wave  function  and  the  final  wavefunction  the  Scatter¬ 
s' 

ing  Matrix  Tfi  can  now  be  obtained. 


v^}  is  the 
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IV.  The  Scattering  Matrix,  Tfi 

This  chapter  develops  and  evaluates  the  scattering  matrix,  T^, 
for  use  in  the  differential  cross  section.  The  total  scattering  matrix 
is  derived  and  each  term  discussed  in  terms  of  particles  exchanged. 

The  spin  states  of  the  initial  and  final  configuration  are  averaged 
and  summed  over,  respectively,  to  obtain  the  scattering  matrix,  Tfi. 
Finally,  the  terms  in  the  scattering  matrix,  T^,  are  evaluated  and 
discussed. 


Scattering  Matrix,  Tfi 

The  total  scattering  matrix,  ,  is  developed  using  the  first 
term  in  the  Born  Series  and  the  initial  and  final  wavefunctions  and 
scattering  potential  developed  in  the  previous  chapter.  The  scattering 
matrix  is  given  by 


Tf.  =  <f I V | i> 


and  using  the  identity, 


<uv.(si>  V(si>  »  =  «v.v! 
1  1  11 


(82) 


(83) 


where  <$v  ,,,  is  the  Kronecker  delta,  the  scattering  matrix  is  given  by 
vi'vi 

Tfi  =  { Bv, ,v-  5i,v£  5-i,v3  -  4v,,v;  s-i,v’2  «i,vd(',2,3|v|1,2,3)  ♦ 

5l,v$  «-!,„•  -  \,V-  «i,v.]<3,1,2|V|1,2,3>  + 

Pv,,»;  si,^  U,vi  -  Vi  V  ^,V.](1.3,2|V|1,2,3)  * 

C!V,,V^  SJ.v$  -  5v,,v^  S-l,v$  <Sl,v'3(2>1’3|V|l,2,3)  + 
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-i.v; 


-  6 


v1’v3  _2,v1 


6^v,](3,2,l|V|l,2,3)} 


»Vi 


M*| 


J2  w2' 


,} 


(84) 


where  (i,j,k  V  1,2,3)  represents  the  spatial  integration  of 
(i,  j,k|v|  1,2,3)  =  y*4'H(  N,i,a)^R(nl,j,B)^e(k)-V- 

4>h(1s,1)'I'h"(1s2,2,3)  dT  (85) 

The  terras  in  the  scattering  matrix  can  be  defined  in  terms  of  ex¬ 
changed  particles.  This  leaves  six  basic  scattering  matrices.  In 
the  case  of  no  exchange,  or  the  direct  scattering  matrix,  T  ,  the 
scattering  matrix  is  given  by 


T ,  = 


[V; 


Jt,v£ 


$1  v'  ”  V*  <^-1  v'  u'  °u  u*  ^ 

s,v3  v1,v1  2,v2  5,v3  p1,p1  P2,y2 


(1,2,3  +  1 ,3,2 |V| 1 ,2,3) 


(86) 


The  scattering  matrix  that  accounts  for  exchange  of  electrons  between 
the  final  H  target  and  H  resultant  atoms  is 


Dv  vi  5_i  v» 

v1,v2  2,v3  2,v1 


"  6^i»V2  6“^,v3  52’vi^  (6 


y1,y1 


6  ,) 

y2»v*3 


(3,1,2  +  2 , I , 3 | V 1 1 ,2,3) 


(87) 


The  scattering  matrix  that  accounts  for  exchange  of  electrons  between 
the  H  target  atom  and  the  free  electron  is 


!,f  =  &vlfU 


'  6-i,V^  'V^V. 


-  6  -  -»  6-£,v’  5i,v 


1  -2»v2  V1 ’  3 
(2,3,1  +  3,2, 1 | V| 1 ,2,3) 


,]*(«. 


y1,y1 


y2,y2 


(88) 


The  scattering  matrix  that  accounts  for  exchange  of  hydrogen  atoms 
between  the  final  H  target  atom  and  H  resultant  atom  is 

th  =  t5Vv-  Ss,v^  6-},vj  -  5vrv^  sii2,n|1- 

(3,1,2  +  2, 1 ,3|v| 1 ,2,3) 


(89) 
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The  scattering  matrix  that  accounts  for  exchange  of  protons  in  the 
hydrogen  atoms  and  electrons  between  the  target  atom  and  free  electron  is 


lp,f 


=  ^vrv3  6^,vi  6-2’v2  "  6vrv3  5_^,v 


5  i 

“2 

(2,3,1  +  3,2, 1 |V| 1 ,2,3) 


+2 


4 


(90) 


The  scattering  matrix  that  accounts  for  exchange  of  protons  between  the 
final  H  atom  is 


=  1"<5  ,5,  ,  -  6 

Lvrv]  ^,v4  ~2 ,  Vq  v 


2  z » v3  vi  »  v2  2,v2 

1,2,3  +  1 , 3 ,2 1 V  j 1,2,3) 


i  «,»  v»3’ 1 

”2»V5  2>v2 


6  ,  6 

V»1 1^2  P 


2’H 


(91) 


In  all  cases  the  last  term  represents  the  spatial  integration  of  the 
wave  function  as  shown  in  Equation  84.  The  total  scattering  matrix 
is  now  written  as 


(92) 


Averaging  and  Summing  Over  Spin  States 

The  scattering  matrix,  Tfi,  depends  explicitly  on  the  spin  of  the 
particles  before  and  after  the  collision.  To  remove  this  dependence 
it  is  necessary  to  sum  over  the  final  spin  states  and  average  over  the 
initial  spin  states.  Table  2  represents  the  possible  spin  configura¬ 
tions  and  the  associated  scattering  matrix  for  each  spin  configuration. 
In  Table  2,  Tft  represents  the  spatial  integration  of  the  scattering 
matrix  TA,  such  as 

Td  =  (3,1,2  +  2,1 ,3|V| 1 ,2,3)  (93) 


Not  all  spin  states  have  been  included,  but  it  can  easily  be  shown  that 
all  other  spin  states  are  just  representations  of  those  in  Table  2. 

The  scattering  matrix  for  the  differential  cross  section  is  then 


Table  2 

Spin  Configurations  and  Associated  Scattering  Matrices 


INITIAL 

FINAL 

1 T  | 

e~ 

£ 

e~ 

p* 

» 

t 

i 

t 

f 

V1 

v2 

^3 

y1 

y2 

V1 

'fe 

^2 

+ 

+ 

- 

+ 

+ 

+ 

+ 

- 

+ 

+ 

!Td- 

Te,H  " 

tp  *  thI 

+ 

+ 

- 

+ 

+ 

+ 

- 

+ 

+ 

+ 

iTd- 

T 

e,f 

tp  *  V1 

+ 

+ 

- 

+ 

+ 

• 

+ 

+ 

+ 

+ 

|Te,H 

‘  Te,f 

*  Tp,f  "  th 

+ 

+ 

- 

+ 

- 

+ 

+ 

- 

+ 

- 

•Td  - 

Te,H^ 

+ 

+ 

- 

+ 

- 

+ 

- 

+ 

+ 

- 

lTd- 

Te,fl 

+ 

- 

+ 

- 

- 

+ 

+ 

+ 

- 

|Te,H 

-v 

+ 

+ 

- 

+ 

- 

+ 

+ 

- 

- 

+ 

V 

thI 

+ 

+ 

- 

+ 

- 

+  ' 

- 

+  ' 

- 

+ 

'V 

TP,rl 

+ 

+ 

- 

+ 

- 

+ 

+ 

- 

+ 

Ith~ 

T  1 
P,f‘ 

lTfil  =  I  <lTd  "  Te,H  -  Tp  +  THl  +  lTd  -  Te,f  -  Tp  +  TP,fl  + 

lTe,H  -  Te,fl+  TP,f  -  THl  +  lTd  *  Te,Hl  +  Td  “  Tp,fi  + 

lTe,H  -  Te,f I  +  lTp  -  THi  +  lTP  -  Tp,fl  +  lTH  -  Tp,fl> 


The  evaluation  of  each  of  the  scattering  matrices  is  presented  in  the 
next  6  sections. 

Evaluation  of  Direct  Scattering  Matrix,  T^ 

The  scattering  matrix  associated  with  no  exchange  of  particles  is 
the  direct  scattering  matrix  and  is  given  by  the  spatial  terms  of 

Td  =  /ri  /%«  ;  1 . 1 )  (4»r  ( nl ,  2 , 2 )  ( 3 )  +  \|>p(nl,3,2)\pe(2)  }• 

V^(1s;1)^(1s212,3)dT 


(95) 


where  \|>h(N;1,1),  <»R(nl;  j,B),  ^(k),  t|>h(1s;1),  >|>H(  1s2;2,3)  and  V  are 
described  by  Equations  77,  78,  80,  74,  75,  and  72,  respectively. 
Written  explicitly  this  is 

_L  f  n* 

Td  =  exp  (5  +  n))F(-n1,1,^/N)F(-n2,1,n/N)  • 

,  _1  <  *  _  _  _  _ 

{(2tr)  2Nr  x  (nl,2)exp(-ikH*RH)exp(-ike*re)  + 

t(2  )-5Nr  X*(nl,3)exp(-iI^*RH)exp(-ike-re)}  • 

{VP  *  V12  *  V13  *  V  -  V3  -  V12)' 


Nh  exp(-r1)exp(-ik_*RQ  +  1(1^  +  ke)*RQ)* 

Nh  exp(-<xr2  -  Br^)  +  exp(-oir3  -  Br2)  exp(ik_'R_)dT  (96) 

* 

where  all  terras  are  defined  in  Chapter  4,  and  x  (nl,i)  is  the  hydrogen 
wave  function  for  the  nl  state.  Denoting 


Nt 


VN’1}  =  n2  exp(-lj  (5+n))F(-n1,l,5/N)F(-n2,1,n/N) 

4»r(  Is,  1 )  =  Nh  expt-r^ ) 

2 

XH-(1s  ;2,3)  =  Nh-  exp(-ar2  -  Br3)  +  exp(-ar3  -  0r2) 


(97) 

(98) 

(99) 


Equation  96  can  be  written  as 
1-15/2 


Td  = 


(2tt) 


/12 


( 1s^;2,3) 


J  ♦JiNjI)  ’  4>h(1s,1)xh  - 
{X*(nl,2)exp(-i,(kH'R’H  +  ke*Fj )  + 
X^nl^expt-idcn'Rn  +  ke-?e) )}  • 


_ I 


exp  (ik_*(R_-  Rq)  +  i(kg  +  1^)  •  RQ) 
VP  +  V12  ♦  V13  -  V  *  V3'  -  V12  * 


(100) 


-18/2 

The  factor  of  (2ir)  is  from  the  constants  in  Equations  97,  98,  and 

99  to  account  for  normalization  of  the  plane  waves  in  the  integral. 

n'*,  Nh-  =  /- and  Nr  =  it-*  if  nl  =  Is  or 
V  2  IT 


1  -1 

Then  NH  =  it  2,  NH 
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>  1 

NR  =  (32H)  if  nl  =  2s.  Making  a  change  of  coordinates  described  in 
Appendix  A,  Equation  99  can  be  written  as 


Td  = 


{2ft) 


-15/2 


/12 


/ 


^(N;1)4'H(1s;1)XH-ns2;2,3)  ' 


(X  (nl’3)exp(-i^*(-  j  r2  +  |J=1T  r3)) 
exp{-ike- (|  ?2  "  5  ?3>  >  +  X*(nl,2)  ' 


exp(-ikH‘(-  3  r3  +  r2))exp(-ike  -  (|  r3  -  ^  r2))}- 

eXp(ffi--(3IW  (?2  *  F3>  -  OTST  ?1>>  ■ 
exp(i(l^  +  •  <2 ('fie')"  **i )  *exp(iK*  (t2  -  7^)  • 


[IV 


T2  “  2  r1  +  3  r2  +  3  r3 


-1 


3  r2  +  1  r3l 


-1 


*T  ~  T  ^ 
1  2 


+ 'T-  +  i  ri  - 


^T1  “  T2  +  2  r1  +  3  r2  "  3  r3 


-1 


2  - 


I  r2  +  3  r3( 


-1 


"  IT1  "  T2  “  2  r1  ‘ 


1 


It,  -^-Iri  4r2-f?3r1  -  | x  -  ?2  +  i  ^  ‘ 


3  r2  +  3  r3| 


-1 


]• 


dt^  dx2  dr1  dr2  dr3 

vrtiere  K  =  k_  -  k^  -  kg.  To  simplify  the  notation,  the  following 
convention  is  used, 

»#  # 


(101) 


VN;1)  =  ^(N;1)exp(+i(kH  +  Ej  •  r,) 


»* 


X  (nl;i,j)  =  x*(nl;i)exp(i  r±  •  k3  -  1^) 

expdFj  •  <^-§Ee)) 

\j»H(1s;1)  =  4>h( Is;  1  )exp(-i  k_  •  r^ 


(102) 


(103) 


(104) 


26 


and 


;2,3)  =x^j-(1s  ;2,3)exp(ik_  *  2 +  *"3 


(r_  +  r_))  (105) 


Using  this  notation  Equation  101  is  written  as 
,-15/2  f  ,# 


Td  = 


(2  u) 


^2 


-J  ^*(N;1)^(1s;1)Xh-(1s2;2,3)  • 

(X  (nl;3,2)  +  x  tnl;2,3)>* 


0 


-  1  r.  -  ’• 

-I 


-1 


1  -  2  - 


T1  "  t2  “  2  r1  +  3  r2  +  3r3'  +  IT1  "  T2  +  2  r1  "3  r2 


»  3r3l  * 


2  - 

"2 
-1 


-1 


1  ~ 


T1  "  T2  +  2  rl  +  3  r2  "  3  r3l  ~  iT1  "  T2  "  2  ri  ~ 


2  -  1- 
3  r2  +  3  r3 


-1 


I  ^1  "  T2  "  2  r1  +  3  r2  "  3  r3l  ”  IT1  "  t2  +  2  r1 


If 

3  2 


*^3l  I’ 

exp(iK*  W2  -  dx2  d?1  d?2  d?3 


(106) 


The  integration  over  t2  space  is  accomplished  using  the  identify 


/ 


=  i|  exp(iK*r) 

|r  -  rl  r 


(107) 


and  Td  becomes 


T,  = 


Att  ( 2n ) 


-15/2  f  ,* 


/12 


/  ^''(NjD^dsiDxjjds*"^^)  • 

{x'*(nl;2,3)  +  X'*(nl;3,2)}  • 

£exp(i£*{-  ^  ri  +  3  F2  +  3  F3)}  +  exp(iK*(^  ?1  - 
§F2  +  I?3})  + 
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r  F 


—  1—  1—  2—  —  1  — 
exp(iK*(^  r1  +  3  r2  “  3  r3^  "  exP(ik*(~  2  r1  ~ 


exp(iK*(-  ^  r1  +  3  r2  ”  f  r3^  "  exP(iK*(3  r-|  " 

3  r2  "  3f  r35)]  * 

dr^  dr^  d?2  dr^  (108) 


Because  of  the  volume  normalization 


T1  IS 


of  plane  waves  the  integration  of 
(2tt)3  (109) 


Performing  this  integration  and  factoring  terms  Equation  108  becomes 
9 


T  -  4tt(2tt)  2  f 

d  "  /12  K2  J 


1* 


(N;1  )vp^(  1  s ;  1 )  • 


{exp(-  2  ^*r1 )  -  exp(^  !(,r1)}dr1  • 

{XH*(nl;2,3)  +  XH*(nl;3,2)}  *  >fo(1s2;2,3)  • 
exp(iK*-~(r2  +  F3))  -  2{exp(iK*(-  3  ^  +  1  ^  * 


dr2  dr 3  (110) 

The  factor  of  2  in  the  second  integral  is  due  to  the  interchangability 
of  r2  and  r3  before  and  after  the  collision. 

The  form  factor  of  the  hydrogen  atom  in  a  ground  state  excited  to 

^  y_ 

its  N  excited  level  is  performed  in  Landau  and  Lifshitz  (Ref  12:578). 


The  second  integral  is  analytically  done  for  simple  initial  wave 
functions.  If  the  final  state  is  Is  then 


/' 


exp(-^r)  exp(il(*r)  •  —  exp(-r)  dr  = 
8  /H  (y+1)/(  (Y+1)2^  K2)2 


(Ilia) 


and  if  the  final  state  is  2s, 


-1 


exp(-yr)  exp(iK*r)  (32tt)  (2-r)  exp(-r/2)dr  = 

( 6 ( y  +  ^)3  -  3(y  +  |)2  +  (2*y  +  2)K2) 


«Y 


+  i)2  +  K2)3 


(111b) 


Using  these  results,  the  integration  of  Equation  110  is  easily  per¬ 
formed,  and 


Tj  =  2(2tt  U-2-  J ^*(N;1){exp(-iK-|  r^  -  exp(iK^  r.,)} 


/S  K2 


i|^(1s;1)  dr  * 


[F(nl,0H-1,^ 


K  + 


k  +  4  k  - 


3( 1+2e)  -  3  e  3(1+e) 


kH)  • 


F(1S’64K  +  3(Wk-+IkH-f  V  + 


F(nl,&<-1,3>  K  +  3 (i+2e)  k-  +  3  ke  "  3(  1+e)  ^ 


F(  Is, a, 3  K  +  3^1+2e)  k_  +  3  ^  "  3  k( 

F(nl,a+1 ,  -  §  K  +  k-  +  I  ke  "  3?  1+e)  kH) 

f(  is,e,j,  k  +  k_  +  i  §  V  - 

F(nl,g+1,  -  §  K  +  Trrglj  I  +  3V  InJeT  V 

F(1s,o,3  K  +  3(  i+2e)  k-  +  3  “  I  ke}  “ 

F(nl,a*-1 ,3  K  +  3^]2e)  k-  +  3  ^  §  V  ’ 

FI  Is,  B,-  3  ^  +  3(1  +2 £ )  k-  +  I^e  “  i(  lie)  “ 


2  T7 
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F(nl,B+1 ,3  K  +  3{1+2e)  k_  +  3  ^  ‘  3  ke) 

F(  1s»a»"  §  ^  +  3(  i+2e)  +  3  "  3(  Wei  V  ’] 


(112) 


where  the  function  F(nl,A,K)  is  defined  as 


F( 1s,A,K)  = 


8/n  (A+l) 
((A+1)2  +  K2)2 


if  nl  =  Is  and 


F(2s,y,K)  = 


/ST  6(y  +  h3  -  3 ( y  +  i)2  +  (2y  +  2_K2 


n  1>2  „2,3 

{  (Y  +  2^  +  K  } 


if  nl  =  2s. 


Exchange  Matrix  for  the  Exchange  of  Electrons  Between  H  Atoms,  T  „ 

_ _  6  ,n 

The  scattering  matrix  associated  with  the  exchange  of  electrons 
between  the  H  atoms  is  given  by 

Te,H  =  ^  K(Ni2.1)%(nli1,2)**(3)-V- 
\|>H-(1s;2,3)\|>H(1s;1)dT  + 


_L  f  * 


3, 1)^R(nl;1,2)^e(2) -V- 


^H-(1s^;2,3)ipH(1s;1)dx 


(113) 


*  *  #  2 
where  iJ>H(N;i,a) ,  vpR  ( nl ;  j ,  3 ) ,  4^00,  \pH~  ( 1  s  ;2,3)  and  ^(Isrl)  are 

defined  in  Equations  77,  78  or  79,  80,  75  and  74,  respectively.  The 
potential  V  is  different  than  Equation  72  on  the  scattering  potential 
in  the  direct  scattering  matrix,  T^.  The  scattering  potential  for 
exchange  of  particles  is  the  electrostatic  potential  between  the  two 
particles.  This  result  is  proved  by  Goldberger  and  Watson  (Ref  9:156-8) 
and  will  not  be  presented  here.  Therefore  the  potential  for  the 


exchange  of  electron  1  and  electron  3  is 

v  =  vi3  =  I*;  ~h\ 


-1 


and  for  the  exchange  of  electron  1  and  electron  2  is 

V  =  V,2  =  |5f,  -  X2f’ 


(114) 


(115) 


Therefore  T  H  can  be  written  as 

0  f  IJ 

4^(  Is;  1  )4<h“(  1s2;2,3)dx  +‘ 


^(N:3’1 


)iPR(nl;1,2)^e(2)-V13- 


( Is ;  1  )\|^-(  1s^;2,3)dx 


(116) 


Because  of  the  indistinguishability  of  the  electrons  the  two  integrals 
in  Equation  116  are  equa 
one  of  the  integrals  or 


in  Equation  116  are  equal,  therefore  T  „  is  just  double  the  value  of 

0  j  *1 


\j^(N;3,1)^-(nl;1 ,2)\{je(2)  *V13’ 
\|^(  Is;  1  )i^-{  1s2;2,3)dx 


(117) 


Using  the  notation  developed  in  the  previous  section,  specifically 

Equations  97,  98  and  99,  T  „  can  be  written  as 

0)0 


le,H 


(2ir) 


-15/2 


/3 


J  ^(N;3)^H(1s;1)xH-(1s2;2,3)x*(nl;1)  * 
exp(i{k_-(R_  -  Rq)  -  ik^  (R^  -  R0)  - 


*e(re  "  V»  • 


V13  * 


(118) 


A  change  of  variables  is  performed  using  the  transformation  matrices 
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developed  in  Appendix  A  for  this  exchange.  Equation  118  is  now 


(  2tt  ) 
SS 


-15/2  f  *  2 

- -  I  vl>H(N;3)4>H(1s;1)xH-(1s^;2,3)x  (nl;1) 


exp(i{k_-(-  1+e  r1  +  r1  +  1+|£  r2  -  ]+|£  r3)  - 
*V(T+e  rl  +  TTe  r1  “  T+e  r3) 


_  _ I  _  _ 

ke‘(r1  +  r2  -  r3)})  * 


_  _  _1 

r1  -  r^  dr^dr.drpdr^dx 


(119) 


The  integration  over  r^  and  T-spaces  can  be  done  analytically  for 

_  _!  _* 

simple  wave  functions.  The  integration  of  the  r^r^  and  r^  can  not  be 
done  analytically  because  the  resultant  wave  function  orbital  coordi¬ 
nates,  ?H,  depends  upon  r 1 ,  ?2  and  r3,  that  is  the  term  exp ( - | -r ^  + 
r1  -  r3|)  appears  in  the  integrand.  Because  of  this  problem,  the 
effect  of  the  exchange  of  electrons  can  not  be  seen  in  the  differential 
cross  section,  because  this  integral  can  not  be  performed  analytically. 
Further  discussions  on  this  matrix  is  presented  later  in  the  Conclu¬ 
sions  and  Recommendations  section. 


Exchange  Matrix  for  the  Exchange  of  Electrons  Between  the  Target  H  Atom 
and  the  Free  Electron,  T  „ 


The  scattering  matrix  associated  with  the  exchange  of  electrons 
between  the  target  H  atom  and  the  free  electron  is  given  by 

Te>f  =  j4(N;2,1)li) R(nl;3,2)vp*(  1 J  *v* 

4^  ( 1  s ;  1 )  4^—  ( 1  s2 ;  2 , 3 )  dir  + 

\|iH(1s;1)ipH-(1s2;2,3)dT  (120) 


*** * 


f 


where  i|>H(N;i,a),  ^R(nl; j,6) ,  i|>e(k),  i^-ds  ;2,3)  and  4>H(  Is;  1 )  are 
defined  in  Equations  77,  78  or  79,  80,  75  and  74,  respectively.  The 
potential  is,  following  the  convention  presented  in  the  previous 
chapter,  or  depending  upon  the  particles  exchanged.  Once 
again  the  electrons  2  and  3  are  interchangable,  therefore 


e,f 


(2ir) 


-15/2 


/3 


/ 


^1N;2)^H(1s;1)XH-(1s2;2,3)X*(nl;3)  ' 


_  _  __T 


exp{i(k_-(^_  -  rq)  -  kH*(RH  -  Rq)  -  ke-(re  -  R^) )}• 
1^  -  X^"1  •  d?  (121) 


where  4^(N;2),  ^(Isjl),  Xo“ds  ;2,3)  and  x  (nl;3)  are  defined  following 
the  notation  in  Equations  97,  98,  99  and  96b,  respectively.  A  change  of 
variables  is  performed  using  the  transformation  matrices  developed  in 
Appendix  A  for  this  exchange.  Equation  121  is  now  written  as 
>-15/2 


e,f 


(  2tt  ) 


/3 


A h(N»2)V1s,1!  XH-(1s2;2,3)x*(nl,3)  • 


__ 


exp{i(k.-(-  ^  r,  .  r,  -  Jifj  r2  *  ^  r3)  - 


7T  1  -*  ~ 


r1  -  r2  + 


e  — 


rd  - 


V'vfe  ri,)) 


1+e  x3' 

|r1  -  r 1 |  dTdr1dr1dr2dr3 


(122) 


The  integration  over  t,  r^  and  r^  space  can  be  done  analytically  for  a 
simple  H~  wave  function.  However  there  is  an  integral  of  the  form 
y*\^(N;2)\l>H(1s,1)exp(i  •  (--^rk_))  ’ 

exp(i  r’  '  (k_  -  ^  (kH  +  Eg) ) )  ‘  -  r’ |  •  dr^r’ 

(123) 


This  can  not  be  done  analytically  and  therefore  the  effect  of  exchange 
of  electrons  can  not  be  included  unless  a  numerical  integration  of  the 


spatial  coordinates  is  done.  A  further  discussion  of  this  scattering 
matrix  is  presented  in  the  Conclusion  and  Recommendation  sections. 


Exchange  Matrix  for  Exchange  of  Protons  Between  the  H  Atoms,  T 


The  scattering  matrix  associated  with  the  exchange  of  protons 
between  the  H  atoms  is  given  by 

;1,2)4»p(nl;2,1)^*(3)-Vp- 

.2. 


n  ctouiiis  g-Lveii 

t  =  —  /  Vn;' 

p  A 1J  u 


—  ( 

A2j  H 


4>„-{1s  ;2,3)^p(1s;1)  dx  + 

# 


* 


1,2}4)R(nl:3,1}^e(2)*Vp* 


iJjh-{1s  ;2,3)\|>H(1s;1)  dx 


( 124) 


The  same  notation  has  been  used  here  that  was  used  in  the  previous 
three  sections.  Performing  the  change  of  variables  and  rewriting  the 
wavefunctions  using  Equations  97,  98  and  99,  Equation  124  can  be  re¬ 
written  as 


T_  = 


(2t r) 


-15/2 


S3 


J il;*(N;1}^H(1s;1)x*(nl;3)xH-(1s2;2,3)  • 


e  — 


exp{i(k_.(-  ri  .  r;  +  ^  r2  +  ^  r3)  - 


1 1t  r1  +  1+e  r1  +  1+e  r3} 

-1 


_ 


ke'(-TTHrt  *r2),>  '  'r1  ‘rll  ' 


dx  dr1  dr2  dr^  dr^ 


(125) 


The  integration  of  x  and  r^  can  be  done  analytically,  but  once  again 
rR  =  |r1  -  r^  +  r^l  and  with  this  term  in  the  resultant  H  atom  wave- 
function  the  integration  must  be  done  numerically. 


Exchange  Matrix  for  Exchange  of  H  Atoms, 

The  scattering  matrix  associated  with  the  exchange  of  H  atoms 
between  the  final  H  atoms  is  given  by 

Td  =  A*( N;2,2)**(nl;1,m*{3)  -{V12  +  Vp}- 

l|'H(ls;1)V(1s2;2’3)dT  + 

~  J  %(N;3,2)i|>*(nl;1,1)^*(2)  *  {  V13  +  Vp  }  • 

t|^H{  "lsjl  )^h~(  1s2;2,3)dix  (126) 

The  same  notation  has  been  used  that  was  previously  used  in  the  pre- 
ceeding  sections.  Performing  the  change  of  variables  and  rewriting 
the  wavefunctions  following  the  format  of  Equations  97,  98  and  99 
Equation  126  becomes 

/o_\-15/2  f  x  ? 

Th  =  “ - J^u( N;3)ipH(1s;1  )x  (nld^-ds4^,3)  * 

exp{i(k_*  (t2  -  T,  +  t?2  .  ?3)  -  r,)  - 

G  1  1  1  — . 

*%’(T1  "  t2  +  2( 1+e)  r1  +  1  r2  "  3( 1+e)  r3} 

V|F2 -Tfe?3H>  '  [lT1  -  ~2  *  5  ?1  * 

1  -  2  _  -1 

3  r2  "  3  ^  + 

1  -  1  -  -  _11 

IT1  "  t2  "  2  r1  +  1  [r2  +  r3}l  J 

d?1  d?2  dr^  dx^  dx2  (127) 

Performing  the  integration  of  x2  and  x^  space,  using  Equations  107  and 
109,  respectively,  gives 
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_  4tt  (  2tt  )  ~9/2  f 
(k_  + 


^(N;3)4/H(1s;1)x*(nl;1)xH-(1s2;2,3)  • 

exp{i(k-'(3nSn’  (?2  ♦  *y  -  snirr  7i>  - 

V^'Clie)  r1  *■  3  r2  "  3H+c)  r3l  " 


ke'lr2  -Tfer3>»  ’ 

{exp(i(k_  +  1^)  *  (2  ri  +  3  r2  “  3  r3^  + 
exp(i(k_  +  T^)  *  (-  ^  ^  ?2  +  ^  F3) )} 


dr^  d^  dr^  ( 128) 

_  —  -2 

This  scattering  matrix  is  proportional  to  (k  +  k^)  which  is  much 
smaller  than  (k  -  -  k  )  which  is  what  the  direct  scattering 

matrix  is  proportional  to.  Therefore  at  the  large  energies  of  the 
incoming  H~  ion,  the  exchange  of  H  atoms  can  be  ignored. 


Exchange  Matrix  for  Exchange  of  Protons  Between  H  Atoms  and  Exchange 
of  Electrons  Between  the  Target  A tan  and  Free  Electron,  t  ~ 

The  scattering  matrix  associated  with  the  exchange  of  protons 

between  the  final  H  atoms  and  the  exchange  of  electrons  between  the 

final  target  H  atom  and  the  free  electron  is  given  by 

Tp,f  =  J %(N;2,2)^-(nl;3,1)4.*n)  '  <V,2  ♦  Vp>  • 

^(1s;1)i|^-(1s2;2,3)  + 

J 4^(N;3,2)xPH-(nl;2,1)^*C1)  *  {  V,3  +  Vp  }  • 

^(1s;1)4^-(1s2;2,3)  dT  .  (129) 

The  same  notation  has  been  used  that  was  used  in  the  previous  five 
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sections.  Performing  the  change  of  variables  and  rewriting  the  wave- 
functions  following  the  format  of  Equations  97,  98  and  99,  Equation 


129  can  be  written  as 
.-15/2 


(2ir) 


p’f  "  S5 


■ 


4<*(N;3)^(1s;1)x*(nl;2)xH-(1s2;2,3)  * 


exp{i(k_-  (--£_r1  +  rH+1TJEr2-Tigr3)- 
V T+e  rH  ■  TTe  r2  +  r3}  “ 


e  — 


1+e  - 


ke’(r1  "  rH  "  1+e  r2  +  r3))} 


{|r1  -  rH|_1  +  |rH  -  r3|“1}dr1  drH  dr2  dr3  dx  (130) 


Again  this  integration  can  not  be  done  analytically.  How  this  affects 
the  total  scattering  matrix,  T^,  is  unknown  because  there  is  no  way 
to  get  a  handle  on  the  size  or  properties  of  this  integral. 


Sunmary 

Of  the  six  basic  scattering  matrices,  only  two  can  be  evaluated 
analytically.  The  scattering  matrix  for  the  exchange  of  H  atoms  was 
shown  to  be  small  compared  to  the  direct  scattering  matrix.  The  re¬ 
maining  four  scattering  matrices  can  not  be  evaluated  analytically  and 
hence  their  contributions  to  the  differential  cross  section  are  unknown. 
Because  of  the  high  initial  energy  of  the  incoming  H“  ion,  it  would  be 
possible  to  assume  that  the  contribution  of  the  remaining  four  scatter¬ 
ing  matrices  to  small.  Following  up  on  this  assumption,  the  total 
scattering  matrix  is  then 

|Tf.|  =  2|Td|  (131) 

In  order  to  prove  the  validity  of  the  assumption,  the  four  remaining 


integrals  must  be  done.  Using  Equation  131  the  differential  and  total 
cross  section  must  be  found  using  numerical  integration  techniques. 

A  discussion  of  the  numerical  integration  technique  is  discussed  in 
Appendix  B. 


38 


I 


V.  Results,  Conclusions  and  Recommendations 


The  purpose  of  this  chapter  is  to  present  the  results  of  the 
computer  program,  discuss  the  results  and  then  recoranend  further  work 
that  should  be  done.  The  computer  program,  Appendix  B,  was  used  to 
calculate  the  differential  cross  sections  for  different  cases.  These 
results  are  discussed  and  specific  conclusions  are  drawn.  Finally 
there  are  some  recommendations  concerning  further  work  in  this  field. 

Results 

The  differential  cross  section  was  calculated  for  many  different 

cases  and  conditions.  The  energy  for  all  the  cross  sections  was  either 

2  MeV  or  200  MeV.  The  factors  varied  were  the  ratio  of  the  electron 

mass  to  the  proton  mass,  the  target  atoms  mass  and  the  final  states. 

Where  possible  the  results  in  this  study  are  compared  to  Genoni  and 

Wright  (Ref  8:4).  The  differential  cross  section  presented  is  really 

not  The  two  are  related  by 

da  .  da 

d6  =  2Tr  sin  6  dfi 

Comparisons  were  made  for  different  incoming  energies.  Figure  2 
and  Figure  3  are  comparisons  of  differential  cross  sections  for  in- 
caning  energies  of  2  MeV  and  200  MeV,  respectively.  The  differential 
cross  sections  calculated  in  this  thesis  are  designated  F(nl),  where 
nl  signifies  the  final  state  of  the  resultant  H  atom.  Those  calculated 
by  Genoni  and  Wright  are  designated  GW(nl).  In  both  cases  there  is  a 
significant  difference  between  the  present  values  and  the  values  from 
Genoni  and  Wright. 

There  are  two  basic  reasons  for  the  differences  in  the  above 
cross  sections.  The  first  difference  is  that  Genoni  and  Wright  include 


AL  CROSS  SECTION 


Figure  2.  Comparison  of  Differential  Cross  Section  for  2  MeV 
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DIFFERENTIAL  CROSS  SECTION 


the  differential  cross  section  associated  with  ionization  of  the 
target  H  atom  and  the  present  value  does  not.  The  second  difference  is 
the  most  important  with  respect  to  the  magnitude  of  the  differential 
cross  section.  Genoni  and  Wright  have  assumed  that  after  the  collision 
the  electron  will  be  in  the  vicinity  of  the  resultant  H  atom  for  an 
appreciable  amount  of  time.  Because  of  this  assumption,  the  Hamilton¬ 
ian  for  the  initial  H”  ion  system  and  the  final  H  -  e"  system  is  the 
same.  They  then  assumed  that  the  final  wavefunction  must  be  orthogonal 
to  the  initial  H~  ion  wavefunction.  They  orthogonalized  according  to 
the  Schmidt  orthogonalization  procedure  or 

|f>  =  \f>  -  <i|f>|i> 

where  |f|  is  the  wave  function  of  the  H  -  e~  system  and  |i>  is  the  H“ 
ion  wave  function.  By  introducing  this  type  of  wave  function  there  is 
a  major  difference  because  of  the  subtraction  of  part  of  the  wave 
function.  This  accounts  for  the  major  differences  in  the  compared  data. 
A  minor,  as  yet  to  be  investigated,  effect  could  possibly  be  that 
Genoni  and  Wright  have  assumed  that  the  ratio  of  the  electron's  mass 
to  the  proton's  mass  is  zero. 

The  effect  of  setting  the  ratio  of  the  electrons  mass  to  the 
protons  mass  to  zero  is  presented  in  Figure  4.  The  ratio  e  was  set 
equal  to  zero  for  the  nl  =  Is  and  2s  and  for  energies  2  MeV  and  200  MeV. 
The  differential  cross  section,  for  e  =  0,  is  designated  F(E,  nl,  0), 
where  E  is  the  incoming  energy  and  nl  is  the  final  state.  When  e  i  0, 
the  cross  section  is  designated  as  F(E,  nl).  It  is  evident  that  if 
e  =  0,  there  is  a  significant  change  of  the  cross  section  at  200  MeV. 
This  is  due  to  the  fact  that  ck_  or  ek^  is  approximately  equal  to  the 
electron  momentum  kg.  Therefore  the  effect  of  the  ratio  e  is  very 


Figure  5.  Differential  Cross  Section  for  D-H  Reaction 
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significant  and  must  be  included  in  the  calculations. 

To  illustrate  the  effect  of  using  heavier  targets,  the  differential 
cross  section  for  nl  =  2s  was  calculated  using  a  target  atom  that  was 
Deuterium.  Because  the  electron’s  orbital  wave  functions  depend  primarily 
upon  the  electrostatic  effects,  the  wavefunctions  for  deuterium  would 
essentially  be  the  same  as  those  of  hydrogen.  Therefore  hydrogen  wave- 
functions  were  used  and  the  mass  of  the  target  was  Mp  =  3673. 1 04  me. 
Figure  5  illustrates  that  the  added  mass  has  decreased  the  size  of  the 
cross  section  but  has  not  significantly  changed  the  shape.  By  in¬ 
creasing  the  mass  of  the  target  atom  the  overall  deflection  of  the  beam 
can  be  decreased,  assuming  a  Hartree-Fock  wave  function  for  the  target 
atom. 

The  total  cross  section  is  calculated  for  the  H  -  H”  reaction. 

Figure  6  compares  the  value  obtained  with  that  obtained  by  Genoni  and 
Wright .  The  total  cross  section  is  high  by  an  order  of  magnitude  all 
cases.  The  reason  for  this  is  because  of  the  larger  differential  cross 
section  obtained  in  this  study. 

Conclusions 

The  results  obtained  in  this  study  compare  in  shape  but  are  con-  • 
sistently  larger  in  magnitude  with  results  obtained  in  earlier  work. 

The  main  reason  for  this  effect  is  due  to  the  difference  in  the  formula¬ 
tion  of  the  scattering  problem.  The  effect  of  exchange  could  not  be 
included  in  the  results  because  of  the  difficulty  in  obtaining  analytic 
solutions  to  the  integrals  involved.  This  effect  could,  in  theory, 
change  the  overall  shape  and  size  of  the  differential  cross  sections 
drastically. 

The  effect  of  heavier  target  atoms  was  investigated  using  a 
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Figure  4.  Effect  of  Setting  me/Mp  =  0 
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Figure  6.  Total  Cross  Section 
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deuterium  atom  instead  of  hydrogen  as  a  target  atom.  The  increase  in 
mass  tended  to  decrease  the  size  of  the  cross  section,  but  the  shape 
was  essential  the  same.  Further  work  with  different  atomic  systems  is 
needed . 

Finally  the  total  cross  section  for  each  final  state  was  larger 
than  that  obtained  by  Genoni  and  Wright.  This  was  due  to  the  larger 
and  wider  differential  cross  sections  calculated  in  this  work. 

In  summary,  it  is  valid  to  use  the  first  term  of  the  Born  Series 
to  calculate  the  differential  cross  section  for  electron  detachment. 

The  incoming  energy  must  be  greater  than  the  binding  energy  or  the 
Born  Approximation  will  need  more  terms  of  the  Born  Series.  Special 
attention  is  called  for  in  regards  to  antisymmetrization  of  the  final 
wave  function,  and  hence  exchange. 

Recommendat ions 

Some  questions  have  been  answered  concerning  the  differential 
cross  section.  However,  there  are  still  sane  areas  that  need  to  be 
investigated.  There  is  a  need  to  understand  how  heavier  targets  affect 
the  differential  cross  section.  This  could  be  investigated  using 
Hartree-Fock  wavefunctions  for  the  target  atom  and  performing  the 
integration  of  the  spatial  coordinates  numerically. 

The  exchange  scattering  matrices  need  to  be  evaluated  and  com¬ 
pared  to  the  direct  scattering  matrix.  By  performing  the  integrations 
numerically  the  assumption  that  the  exchange  integrals  are  small,  com¬ 
pared  to  the  direct  scattering  matrix  can  be  investigated.  This  inte¬ 
gration  would  most  probably  be  done  numerically  using  a  Monte-Carlo 
integration  routine. 
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There  should  be  an  investigation  into  the  effect  of  using  better 
wavefunctions  for  the  H~  ion.  Different  wavef unctions  could  be  used 
to  see  how  the  differential  cross  section  changes.  Once  again  a 
Monte-Carlo  integration  routine  would  be  needed  because  this  calcula¬ 
tion  would  require  a  15-fold  integration.  By  combining  the  above 
recommendations,  a  general  program  could  be  produced  that  would  use  a 
H"  wavefunction  that  has  been  shown  to  be  reasonably  accurate,  and 
given  an  arbitrary  target  atom,  the  program  could  calculate  a  differ¬ 
ential  cross  section  for  this  system  with  exchange  included. 

A  final  recommendation  would  be  to  see  if  it  is  possible  to  sum  the 
Born  Series  and  obtain  a  better  approximation  to  the  scattering  poten¬ 
tial.  This  is  most  probably  not  a  simple  task  and  should  only  be  done 
if  absolutely  necessary.  An  example  of  when  this  would  be  necessary 
would  be  when  the  incoming  energy  of  the  H"  ion  is  very  small  or  on  the 
order  of  the  binding  energy  of  the  H  system. 
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Apf?  .  A 


Development  of  Integration  Variables 

It  is  necessary  to  develop  a  coordinate  system  that  spans  the 
initial  configuration  space,  the  final  configuration  space  and  the 
scattering  potential  space.  A  basis  set  of  coordinates  are  the 
position  vectors  of  each  particle  relative  to  the  origin,  this  basis  is 
designated 

1  =  (%v  Y2,  X3,  R3)  { A— 1 ) 

where  X.  is  the  position  vector  of  electron  i  and  R .  is  the  position 
1  j 

vector  of  proton  j.  The  general  procedure  will  be  done  for  the  case 
of  no  exchange  of  particles,  but  is  easily  extended  to  the  cases  of 
exchange . 

The  initial  configuration  consists  of  electron  1  and  proton  1 
in  the  target  H  atom  and  electron  2  and  3  and  proton  2  in  the  H-  ion. 
This  configuration  is  described  completely  by  the  coordinates 

\  =  (r*1 ,  Rq,  r2,  r3,  R_)  (A-2) 

I’V'i 

where  r^  is  the  position  vector  of  the  i  electron  relative  to  the 
proton  it  is  bound  to,  and  7?0  and  are  the  position  vectors  of  the 
center  of  mass  of  the  H  target  atom  and  H“  ion  respectively.  7 
is  written  in  terms  of  Z 
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Obviously  r1  =  r1  and  rH  =  r^,  therefore  the  above  matrix  has  arank  of 
3,  and  two  independent  vectors  can  be  chosen  that  can  be  used  to 
describe  the  system.  t1  and  are  chosen  so  that  they  are  orthogonal 
to  ,  r2  and  r^  and  normalized  so  that  the  determinant  of  the  matrix 
is  1,  this  gives 
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Now  Z  can  be  written  in  terms  of  Z  or 
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and  writing  Z  and  Z  in  terms  of  2  ,  give 
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•  Z  ( A— 10) 


It  is  now  possible  to  analytically  integrate  over  the  Z*  set. 

The  same  process  was  used  to  develop  an  integration  basis  set 
when  particles  were  exchanged.  The  initial  configuration  in  the  Z 

basis  set  is  the  same  in  all  cases.  The  final  configuration  is  pre- 

_  _? 

sented  in  the  Z  basis  set  and  the  exchange  explained.  The  Z  basis 

set  is  then  presented  and  finally  ZQ  and  ZQ  are  written  in  terms  of  the 
Z  coordinate  set  for  each  exchange. 

Exchange  of  Electrons  Between  H  Atoms 

t  = 

o 

In  this  exchange  electron  1  and  proton  2  are  in  the  resultant  H  atom 
while  electron  3  and  proton  1  are  in  the  H  target  atom  with  electron  2 
free.  The  transformation  matrix  is 
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In  this  exchange  electron  2  and  electron  1  have  exchanged  positions, 
that  is  electron  1  is  free,  electron  2  is  in  the  target  atom  and  electron 


3  is  in  the  resultant  H  atom.  The  transformation  matrix  is 
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Exchange  of  Protons  Between  H  Atoms 

The  final  configuration 

is  given  by 
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(A-19) 


In  this  exchange,  proton  1  is  in  the  resultant  H  atom,  proton  2  is 
in  the  H  target  atom  and  electron  2  is  free.  There  has  been  no  exchange 
between  electron  1  and  3.  The  transformation  matrix  is  given  by 


(A-20) 


where  Z  =  (r^  r\  r2,  r3,  t).  Writing  ZQ  and  in  terms  of  1  gives 
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Exchange  of  H  Atoms 

The  final  configuration  is  given  by 
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In  this  exchange,  proton  2  and  electron  3  are  in  the  H  target  atom, 
proton  1  and  electron  1  are  in  the  resultant  H  atom  and  electron  2  is 
free.  The  transformation  matrix  is 
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where  z  s  (r^,  t  ,  r £,  t^).  Writing  and  1Q  in  terms  of  Z 

gives 
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Exchange  of  Protons  Between  H  Atoms  and  Exchange  of  Electrons  Between 
Target  H  Atom  and  Free  Electron 


The  final  configuration  is 

given  by 
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In  this  exchange,  proton  2  and  electron  3  are  in  the  target  H  atom, 
proton  1  and  electron  2  are  in  the  resultant  H  atom,  and  electron  1 
is  free.  The  transformation  matrix  is  then 
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Appendix  B 


Program  to  Numerically  Integrate  the 
Differential  and  Total  Cross  Sections 


Because  of  the  complexity  of  the  terms  in  the  differential  and 
total  cross  section,  the  integration  over  the  electron  momentum,  kg, 
can  not  be  done  analytically.  A  program  was  written  to  numerically 
integrate  over  the  electron  momentum,  using  a  Gauss-Legendre  integra¬ 
tion  routine.  The  magnitude  of  the  electron  momentum  was  integrated 
over  using  a  50-point  Gauss-Legendre,  while  the  integration  over  the 
electron  angle  was  done  using  a  24-point  Gauss-Legendre  integration 
routine. 

The  incoming  energy  in  the  lab  frame  for  the  H~  ion,  the  final 
states  of  the  H  atoms  and  the  scattering  angles  are  inputted  to  the 
program.  The  program  then  transforms  the  configuration  to  the  center 
of  mass  frame  via  relativistic  mechanics.  The  Gauss-Legendre  integra¬ 
tion  routine  integrates  over  electron  momentum  and  angle  and  then 
transforms  the  differential  cross  section  to  the  laboratory  frame  using 


where  r  =  V  /V~,  V„  is  the  velocity  of  the  final  resultant  H  atom, 
cm  i  t 

The  total  cross  section  is  then  found  by  using  a  trapezoid  integration 
over  the  differential  cross  section. 

The  integration  routine  was  tested  by  comparing  the  results  of  the 
Gauss-Legendre  integration  to  the  results  computed  via  Rhomberg  integra¬ 
tion  routine.  The  results  of  both  integrations  agreed  to  1  percent. 

This  signifies  that  either  both  of  the  integration  routines  are  wrong 
or  both  are  right.  The  probability  that  both  are  wrong  but  give  the 


same  result  is  very  small,  so  it  can  oe  assumed  that  the  integration 
routine  is  correct. 

There  is  no  simple  analytic  solution  that  can  be  easily  tested 
using  this  program.  However,  each  subroutine  and  function  routines 
were  tested  with  analytic  solutions  where  possible  to  insure  that  each 
was  working  correctly. 

A  listing  of  the  program  follows.  Program  MAIN  contains  all  the 
kinematic  transformations  and  the  integration  routine.  Subroutine  FION 
consists  of  the  integral  of  the  differential  cross  section.  Function 
FF  calculates  the  atomic  form  factor  for  a  general  excited  state  N  of  a 
hydrogen  atom.  Function  F  calculates  the  atonic  form  factor  for  the 
hydrogen  atom  transitioning  to  the  Is  or  2s  states. 
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